function xi=construct_xi_ext(x,n,k,factor_extend)

a = min(x)-0.001;
b = max(x)+0.001;

xi_n = linspace(a,b,n+1)';

xi_neg = xi_n(1,1) + factor_extend*(b-a)*(-k:1:-1)';
xi_pos = xi_n(n+1,1) + factor_extend*(b-a)*(1:1:k)';

xi = [xi_neg; xi_n; xi_pos];

    %
